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^ : Abstract 

li I This is a concise review of S. Chandrasekhar's research contributions to astrophysics, ranging 

[-^ ■ from his early studies on white dwarfs using relativistic quantum statistics to topics as diverse as dy- 



namical friction, negative hydrogen ion, fluid dynamical instabilities, black holes and gravitational 
waves. The exposition is based on simple physical explanations in the context of observational as- 
tronomy, addressed primarily to the undergraduate students. Black holes and their role as central 
engines of active, compact, high energy sources have been discussed in some details. 



I INTRODUCTION 

The impactful research journey of Subrahmanyan Chandrasekhar began on July 31, 1930, 
from Bombay port on a ship. The 19 year old Chandra was on his way to England for higher 
studies. Armed with his understanding of Fowler's work on white dwarfs^, Chandra was 
immersed in the mathematical equations describing these dense objects, during that voyage. 
He had realized that Fowler's theory needed modification, since for sufficiently massive white 
dwarfs, particle number densities could be so high that a large fraction of electrons would 
be occupying very high energy levels, moving with relativistic velocities. 

At this juncture, a quick summary of stellar evolution theory is in store. In main sequence 
stars (like Sun), nuclear fusion of hydrogen to helium supplies the required thermal energy to 
stall gravitational contraction of a star, enabling it to attain a quasi-hydrostatic equilibrium. 
As the star advances in age, a further sequence of nuclear fusion reactions gets activated in 
its core - helium burning to carbon and oxygen, carbon burning to sodium and magnesium 
and so on, if the star is massive enough, till the formation of iron-rich core. Iron nucleus 
being the most stable one, further nuclear burning is energetically not feasible. As the 
iron core cools, the thermal pressure due to the left over heat energy is insufficient to halt 
the collapse under its own weight. Eventually the electron density becomes so high that 
the electron degeneracy pressure comes into play, and prevents further contraction. For 
low mass stars however, white dwarf stage is reached with the formation of a degenerate 
carbon-oxygen core after the red giant phase resulting from helium burning. 

Degeneracy pressure is a consequence of quantum statistics in extremely dense matter. 
Pauli exclusion principle (PEP) states that no two identical fermions can have the same 
state. Electrons, protons, neutrons, neutrinos, etc., being spin half particles, are fermions. 
According to PEP, in a gravitationally bound system like the iron-rich core of an evolved star, 
all the electrons cannot occupy the lowest energy level (unlike, what happens to identical 
bosons in Bose-Einstein condensates, e.g. He-4 superfluid). So, the energy levels are filled 
up with two electrons (one with spin up state and the other with spin down) per orbital, as 
demanded by the PEP. More the density of electrons, higher is the electronic energy level 
that gets to be occupied. 

Gravitational shrinking of such a dense core leads to an increase in electron density, 
thereby encountering a stronger resistance since the contraction implies putting electrons 



at higher energy levels. Therefore, in such a 'degenerate matter' dominated system, the 
gravitational collapse instead of lowering the total energy of the star tends to increase it. 
The ensuing pressure against shrinking, arising out of PEP in such an electron-rich dense 
matter is called electron degeneracy pressure (EDP). A white dwarf is a star that is in 
hydrostatic equilibrium not because of thermal pressure but due to the EDP that counteracts 
gravitational contraction. Fowler had assumed electrons to be moving non-relativistically 
inside the degenerate core, and had arrived at the result that the EDP of a white dwarf is 
proportional to p^^^, where p is the density of the core^. 

II CHANDRASEKHAR LIMIT AND COMPACT OBJECTS 

Having realized that electrons would be randomly whizzing around in the enormously 
dense core with relativistic speeds, Chandra incorporated Einstein's special relativity in the 
analysis of white dwarfs, and found that the EDP is actually proportional to p^'^ instead, 
demonstrating that the relativistic degeneracy pressure does not increase as rapidly as in 
Fowler's case. Employing a physically realistic analysis of the relativistic problem of a dense 
star ruled by a polytropic equation of state, in which gravity is countered by the EDP, he 
arrived at the celebrated Chandrasekhar mass limit^, 

0.2 /;ic\3/2 

where %, G, c, rup and pe are the reduced Planck's constant, Newton's gravitational constant, 
speed of light, mass of a proton and mean molecular weight per electron, respectively. It is 
remarkable that such a significant result concerning stars should be expressible in terms of 
fundamental quantities (except for /ig). In white dwarfs, the value of /ig is about 2, so that 
from eq.(l) one finds the limit to be Mch ~ 1-4 Mq, where Mq = 2 x 10^*^ kg is the Sun's 
mass. 

During the voyage, Chandra was unaware that Anderson in 1929 and Stoner in 1930 had 
independently applied special relativity to obtain mass limits for a degenerate, dense star of 
uniform density but without taking into account the condition of hydrostatic equilibrium^'^'^. 
When Chandra showed his work to Fowler after reaching Cambridge, Fowler pointed out 
contributions of Anderson and Stoner to him. Chandra added these references to his papers 
on relativistic degeneracy in white dwarf stars ^. Incidentally, Landau too had arrived at a 



mass limit independently in 1931, which appeared in print one year later^. 

Chandrasekhar mass limit entails that no white dwarf with mass greater than this limit 
can hold out against gravitational collapse. All the white dwarfs discovered so far (e.g. Sirius 
B, the companion star to Sirius), have mass less than Mch- For niasses beyond this limit, 
two prescient ideas were put forward independently, that played important roles later - one 
of Landau^, before the discovery of neutrons by Chadwick in 1932 and the other by Baade 
and Zwicky®'^, immediately after the discovery. Landau had speculated that for stellar cores 
whose mass exceeded Mch-, the density would become so large due to shrinking that the 
atomic nuclei in the core would come in contact with each other - the whole core turning 
into a giant nucleus^. Baade and Zwicky, while attributing the origin of cosmic rays to 
stellar explosions called supernovae, correctly identified the energy liberated due to sudden 
decrease in the gravitational potential energy (as the core collapses rapidly to form a neutron 
star of radius ~ 10 km) as the one that powers supernova explosion^'^. A core with mass M^, 
shrinking from a large size to a neutron-rich ball with radius Re-, has to give up an energy, 

since its gravitational potential energy decreases to ~ —E^xp- For a 1.4 Mq core collapsing 
to form a neutron star of radius i?c ~ 10 km, the energy E^xp available for explosion is as 
high as ~ 10^^ ergs. 

Why does the core become neutron-rich? As the core shrinks, its density rises till it 
reaches nucleonic values ~ 10^^ - 10^^ gm/cm^, when protons in the core transform into 
neutrons by capturing electrons and emitting neutrinos^°. Neutrinos, being weakly interact- 
ing particles, escape from the core. While in the neutron-rich core, the neutron degeneracy 
pressure (arising from PEP, as neutrons too are spin half particles) prevents further gravi- 
tational contraction, resulting in the formation of a neutron star. 

With the detection of periodic emission of radio-pulses from a source by Jocelyn Bell and 
Anthony Hewish in 1967, existence of neutron stars as pulsars was established. Pulsars are 
rapidly spinning neutron stars with rotation period ranging from about few milli-seconds 
to few seconds. The observed pulses are due to electromagnetic radiation from accelerated 
charge particles moving along strong magnetic field lines inclined with respect to the rotation 
axis (The polar magnetic field strengths vary from ~ 10^° to ~ 10^"^ gauss). Recently, 
a milli-second pulsar was found to have a mass of ~ 2 Mq, determined using a general 



relativistic effect called Shapiro delay in which radiation grazing past a compact, massive 
object, arrives at the observer with a time lag because of the strongly curved space-time 
geometry it encounters near the massive star^^. 

As long as the core is lighter than about 2 — 3 Mq, it can survive as a neutron star 
(The mass limit in this case is uncertain as it depends crucially on the equation of state for 
nuclear matter which, for such huge densities existing inside neutron stars, is unknown^^'^^). 
The released neutrinos, after travelling long distances, eventually lose their energy to the 
stellar envelope, causing the latter to be blown apart, giving rise to a Type II supernova. 
Measurements concerning detected neutrinos from the supernova SN 1987A indicate that 
these ultra-light, weakly interacting particles carry away 99% of the gravitational binding 
energy released from the collapsing core, lending credence to the neutrino driven explosion 
models^'^. 

The observed masses for neutron stars do not appear to exceed ~ 3 Mq^^'^'^, suggesting 
that a massive star whose core is heaver than this limit, would certainly collapse to form 
a black hole. The long duration gamma ray burst sources that exhibit prompt gamma 
emissions with photons having energy predominantly in 0.1 - 1 MeV range, and lasting for 
about 2 - 1000 s are likely to be massive cores collapsing to form black holes^^. Eddington 
was extremely critical of Chandra's work which suggested that a star exceeding the mass 
limit would shrink gravitationally to a point singularity^^. Three decades later, Penrose and 
Hawking, using Raychaudhuri equation, proved the seminal singularity theorems. According 
to these theorems, gravitational collapse of normal matter generically lead to formation of 
singularities, namely, creation of black holes^^"^^. 

Ill DYNAMICAL FRICTION 

Chandra played a stellar role in the research area of gravitational dynamics of point 
masses, from 1939 to 1944, that culminated in his celebrated papers on dynamical 
friction^^'^^. There is abundance of gravitationally bound systems of massive objects in 
the cosmos like globular clusters, galaxies, clusters of galaxies, etc. Entities that make up 
these bound systems, apart from moving in self-created gravitational potential wells, also 
suffer two-body gravitational encounters, resulting in exchange of energy and momentum. 
Chandra discovered that a massive body in motion, surrounded by a swarm of other less 



massive objects, suffers deceleration that is proportional to its mass^^. 

This dynamical friction originates from cumulative gravitational encounters that the mas- 
sive body experiences due to the presence of other objects in the background. The physical 
cause of dynamical friction can be intuitively understood by going to the reference frame in 
which the more massive body is at rest. In this reference frame, the swarm of background 
objects moving past the massive body get focussed behind it because of its attractive force, 
forming a wake of higher number density of objects, and hence of higher mass density. 
When we switch back to the original frame in which the massive body is moving, we find 
that the mass density of the wake behind the larger mass is greater than the density ahead. 
Consequently, it experiences a greater gravitational pull from behind, and thereby suffers 
a gravitational drag force whose magnitude is proportional to the square of its mass and 
inversely proportional to the square of its speed^°'^^. 

Observational support for dynamical friction comes from sinking of globular clusters 
towards the central regions of galaxies, and galactic cannibalism in which the orbit of a 
satellite galaxy decays, leading eventually to its merger with the bigger galaxy^^'^^. 

IV NEGATIVE HYDROGEN ION 

During the early forties, Chandra was also involved with the quantum theory of negative 
hydrogen ion. Now, can a proton capture two electrons to form a charged bound state? How 
is it relevant to astrophysics? The first issue had been settled by Bethe in 1929 who showed 
that quantum mechanics indeed predicts formation of H~ ions^^. As to the second question, 
it has been found over the years that H^ is a weakly bound system with a binding energy 
of ~ 0.75 eV. Since it takes only about 0.75 eV to knock off the extra electron from H~ , 
its life-time under terrestrial conditions is small but in thin and tenuous plasma where the 
collision frequency is low, one expects negative hydrogen ions to survive for longer duration. 

Early on, Wildt had foreseen that H~ would form in large numbers in the upper atmo- 
sphere of Sun, because of the abundant presence of hydrogen atoms and electrons there. He 
had also come to the conclusion that photo-detachment of H' would contribute greatly to 
solar opacity, since radiation from Sun would be attenuated as they photo-ionize H^ ions 
on their way out^"^~^^. 

Chandra and his collaborators played a significant role in calculating H~ photo- 
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absorption matrix element, so crucial for estimating the quantum probability (and, therefore, 
the cross-section) of photo-ionization of negative hydrogen ion^^^^^. The opacity or the op- 
tical depth is proportional to the photo-absorption cross-section a as well as n, the number 
density of iJ~. This is because, the number of photo-ionizations per photon per unit time 
is c n (7, so that the mean free path length for photons is simply, 

n a 

The optical depth essentially is the ratio of the geometrical path length traversed by the 
radiation to mean free path length / (i.e., it is the number of absorptions suffered by the 
photons on an average). 

The negative hydrogen ion has only the ground state as a bound state. There is no 
singly excited state that is a bound state. As a result, photons with energy above 0.75 eV, 
executing random walks out of Sun due to multiple scatterings, would be absorbed by H~ 
ions after detaching their extra electrons to the continuum. This is the major cause for solar 
opacity in the infra-red to visible range of the electromagnetic spectrum. 

In 1943, Chandrasekhar and Krogdahl drew attention to the fact that the dominant 
contribution to the matrix element describing the probability amplitude for such a photo- 
dissociation comes from the wavefunction at large distances (several times the Bohr radius). 
Therefore, an accurate knowledge of electronic wavefunction of H^ is required to compute 
the matrix element correctly^^. 

Chandra and his collaborators made seminal contributions in calculating the continuous 
absorption coefficient kx as a function of the photon wavelength A, taking into account 
dipole-length and dipole-velocity formulae, that provided a solid theoretical foundation for 
the characteristic k\ - \ plot which exhibits a rise in the range 4000 to 9000 angstroms and 
then drops to a minimum at 16000 angstroms, with a subsequent rise^^. 

The negatively charged hydrogen ion also plays an important role in cyclotrons and 
particle accelerators^^. The advantages in making use of H~ arise out of the possibility of 
accelerating them by applying electric fields and obtaining hot neutral beams in Tokamaks 
(like in ITER)^^. This is because of the relative ease in detaching its extra electron when 
H" ion is present in the gas cells. 



V MAGNETOHYDRODYNAMICS 

Astrophysical entities are usually permeated with magnetic fields, be it planets like earth, 
Jupiter, etc.. Sun, sunspots, stars, flares, spiral arms of Milky Way, galaxies, and so on. Mag- 
netic field in a conducting medium like metal or plasma wears out due to Ohmic dissipation. 
How does terrestrial magnetic field, presumably generated by the electric currents flowing 
in the molten, conducting and rotating core of Earth, counter Ohmic decay? 

Dynamo theories involving differential rotation and convection in conducting fluids are 
invoked to solve this riddle. However, Cowling had proved that magnetohydrodynamical 
flows with axisymmetric geometry will always entail a decaying magnetic field'^^. About 
two decades later. Backus and Chandra generalized Cowling's theorem^^. In this context, 
Chandra studied the possibility of lengthening the decay duration so that an axisymmetric 
dynamo provides a feasible explanation for geomagnetism^^. It was immediately followed by 
a paper in which Backus showed that the increase was not large enough to be of geophysical 
interest^". Chandra studied several fluid dynamical stability problems employing variational 
methods that have interesting consequences'^^ '*^^. 

An evolved binary system, consisting of a Roche lobe^*^ fllling star, spewing out gaseous 
matter, and a massive compact object (MCO) like a neutron star or a black hole (BH), both 
going around the common centre of mass, very often acts as a luminous source of high energy 
photons. In such a binary system, gas leaking out from the bloated star cannot radially fall 
on the MCO as it has angular momentum. Instead, it spirals inwards, forming an accretion 
disc around the MCO so that each tiny gaseous volume element of the disc moves along a 
circular Keplerian orb it '^^. 

For a thin disc with a total mass much less than the mass M of the MCO, the Keplerian 
speed v{r) of a fluid element at a distance r is given by. 



v{r) = ^^- , (3) 

Eq.(3) implies that the fluid elements of the accretion disc rotate differentially. Farther the 
element from the MCO, lower is its circular speed. Differential rotation leads to viscous 
rubbing of neighbouring fluid elements at varying distances, causing the accretion disc to 
heat up. If the disc is sufficiently hot, it emits copious amount of electromagnetic radiation 
with a spectrum ranging from visible wavelengths to UV photons and X-rays. 



There are strong observational evidences that the rapidly time varying, intense X-ray 
sources, like Cygnus X-1, are accreting black holes (see section VII). Essentially, the gravi- 
tational potential energy of the gas spiralling in, gets converted into radiative energy at the 
rate corresponding to a luminosity of, 

GMrh 



min 



(4) 



where m, rmm and e are the rate of mass accretion, minimum distance reached by the 
infalling gas and an efficiency factor for the conversion of gravitational energy to radiation, 
respectively. The importance of accretion on to compact objects is evident from eq.(4), since 
source luminosity is larger for smaller values of Tmin- Similarly, a luminous source requires 
larger rates of accretion and higher conversion efficiencies. 

For the efficiency factor e to be large, the accretion disc is required to have a high viscosity. 
The physics of the mechanism responsible for large viscosities in the disc is an active area of 
research. Interestingly, as shown by Balbus and Hawley in 1991, the Chandrasekhar insta- 
bility might be the key to the origin of accretion disc viscosity'^^ . Chandra had pointed out 
that a differentially rotating, conducting and magnetized incompressible fluid in a cylindrical 
configuration, is unstable with respect to oscillating axisymmetric perturbations^^. 

While investigating Rayleigh-Benard convection in conducting and viscous fluids threaded 
with magnetic fleld, Chandra studied the onset of convection and its dependence on a di- 
mensionless number Q, representing the square of the ratio of magnetic force to viscous 
force*^^. Today, this number Q is referred to as Chandrasekhar number (or, also as the 
square of Hartmann number). Chandra made several other contributions in the fleld of 
plasma physics and magnetohydrodynamics that had far reaching consequences^®. 

VI CHANDRASEKHAR-FRIEDMAN-SCHUTZ INSTABILITY 

While studying self-gravitating and rotating fluid conflgurations, Chandra showed that 
a uniformly dense and uniformly rotating incompressible spheroid is unstable because of 
non-radial perturbations, causing emission of gravitational radiation*^^. According to Ein- 
stein's general relativity, the curvature of space-time geometry manifests as gravitational 
force. Gravitational radiations are wave-like perturbations in the space-time geometry that 
propagate with speed of light, general relativity being a relativistic theory of gravitation. 
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Gravitional waves are radiated whenever the quadrupole moment of the mass distribution 
in a source changes with time. Friedman and Schutz extended Chandra's findings in 1978, 
and demonstrated the existence of gravitational wave driven instabihty in the general case 
of rotating and self-gravitating stars made of perfect fiuid^'^. 

An intuitive way to comprehend this Chandrasekhar-Friedman-Schutz (CFS) instability 
is to look at a particular perturbation mode in a rotating star that is retrograde, i.e. moving 
in the backward sense relative to the fiuid element going around. According to general 
relativity, the space-time geometry around a rotating body is such that inertial frames are 
dragged along the direction of rotation (This has been recently verified by the Gravity Probe 
B satellite-borne experiments^). The frame dragging, therefore, would make the retrograde 
mode appear prograde to an inertial observer far away from the star. Gravitational waves 
emitted by this mode will carry positive angular momentum (i.e. having the same sense as 
the angular momentum of the fiuid element) as measured in the distant inertial frame. Since, 
the total angular momentum is conserved, gravitational radiation carrying away positive 
angular momentum from the mode, would make the retrograde mode go around more rapidly 
in the opposite direction, leading to an instability. 

In 1998, Andersson demonstrated that a class of toroidal perturbations (the so called 
r-modes) in a rotating star are generically unstable because of the gravitational wave driven 
CFS instability s^. It was immediately followed by papers arguing that the r-mode instability 
would put brakes on the rotation of a newly born and rapidly spinning neutron star^*^'^^. 
Consequently, as the neutron star spins down, a substantial amount of its rotational energy 
is radiated away as gravitational waves, making it a likely candidate for future detection 
by the laser interferometric gravitational wave detectors, namely, the LIGOs^^'^^. The CFS 
instability may soon be put to experimental tests. 

VII BLACK HOLES AND GRAVITATIONAL WAVES 

Chandra called the astrophysical BHs the most perfect macroscopic objects^^. Macro- 
scopic entities - like chairs, books, computers, etc. around us, need an astronomically large 
number of characteristics each for their description. A sugar cube, for instance, would need 
not only its mass, density, temperature, but also amount and nature of trace compounds 
present, the manner in which sugar molecules are stacked, porosity, surface granularities. 
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etc, for its proper specification. In contrast, a BH is cliaracterized by just tliree pliysical 
quantities - its mass, cliarge and angular momentum. 

Scliwarzscliild BHs do not possess charge or angular momentum, while Kerr BHs rotate 
but have no charge. On the other hand, Reissner-Nordstrom BHs have charge but do not 
rotate. Kerr-Newman BHs are theoretically the most general ones, as they possess non-zero 
mass, charge and angular momentum. Astrophysical black holes are all likely to be Kerr BHs 
since charge of a BH would get neutralized by the capture of oppositely charged particles 
present in the cosmic rays and other ambient matter, and since most cosmic objects possess 
angular momentum. Chandra was particularly fascinated by the stationary, axisymmetric 
vacuum solutions of Einstein equations that described the Kerr BHs. 

BHs are characterized by a fictitious spherical surface called the event horizon centred 
around the point singularity created by the collapse of matter. Nothing can escape from 
regions enclosed within the event horizon. For a Schwarzschild BH of mass M, the radius 
of the event horizon is given by the Schwarzschild radius, 

2GM fi/ M \ , 

But do BHs exist in the real universe? Classical BHs by themselves do not radiate. 
Hawking radiation, which is of quantum mechanical origin, from astrophysical BHs, is too 
miniscule in amount to be of any observational significance^^. So, how does one find BHs in 
nature? In conventional astronomy, their detection relies on the presence of gas or stars in 
their vicinity and the ensuing stellar or dissipative gas dynamics around an accreting MCO. 
As discussed in section V, if the MCO has an accretion disc around it like in galactic X-ray 
sources, quasars, blazars or radio-galaxies, the swirling and inward spiralling gas gets heated 
up, emitting radio, optical, UV and X-ray photons, often accompanied by large scale jets^®. 

If gas can spiral down to a distance Tmin = oi Rs from the central BH, then according to 

eqs. (4) and (5) the radiation luminosity is given by, 

0.5e . 2 /^\ 

L = mc . (6) 

a 

The real parameter a quantifies the proximity to the central BH. Eq.(6) tells us that accretion 
taking place close to the event-horizon can convert an appreciable fraction of rest energy mc^ 
of the infiowing gas. Higher the accretion rate m, larger is the luminosity L. (Provided that 
fiuid viscosities in the disc are large enough to give rise to higher efficiencies e, as discussed 
in section V.) 
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The central engine for a quasar or a blazar is, in all likelihood, an accreting supermassive 
BH with M lying in the range 10^- 10^ Mq^^ . Invoking eq.(6) with sufficiently large accretion 
rates for blazars, one can theoretically explain high luminosities (at times, exceeding 10''* 
erg/s) observed in these sources. 

Quasars and blazars also exhibit fluctuating X-ray luminosities on time scales of only 
few hours. One can derive an upper limit for the size of the central engine from causality 
arguments. If the observed time scale over which the luminosity varies accreciably is At, 
the size of the source participating in emission of photons cannot be larger than cAt. This 
is because, ffistly, every part of the entire region must be causally connected to each other 
and, secondly, special relativity tells us that parts of the region can physically communicate 
with each other (to remain in causal touch) only with speeds < c. X-ray variability on time 
scales of an hour corresponds to a causal size < 10^ km. Now, from eq.(5), a BH of mass 
3 X 10* Mq has a Schwarzschild radius of about 10^ km. Short time fluctuations and central 
engines involving gas dynamics close to the event horizon of BHs, fit together neatly. 

Observational evidence for accreting super-massive BHs comes not only from short time 
variability of X-ray fiuxes but also from the details of the continuum spectra (e.g. presence 
of the big blue bump in quasar spectra) observed in these active sources. Hence, quasars, 
blazars and powerful radio-galaxies are most probably distant galaxies housing acccreting 
supermassive BHs with mass in excess of 10^ Mq in their central regions^^. 

Similarly, by monitoring stellar dynamics around the central region of Milky Way for 
decades, one infers that the Galactic nucleus contains a heavy and compact object, most 
likely to be a supermassive BH with a mass of about 4 x 10^ Mq, within a radius of 10^^ 
km from the Galactic Centre^^. It is interesting to note that the Chandra X-ray observatory 
(launched on July 23, 1999, and named after S. Chandrasekhar) revealed the presence of a 
X-ray source as well as hot gas with high pressure and strong magnetic field in the vicinity 
of the Galactic Centre. 

However, these are indirect detections, implying strictly speaking the presence of a very 
compact, massive central object. Inference of an astrophysical BH, although very likely, 
relies on theoretical interpretation. What happens when a BH is perturbed by incident 
gravitational waves or electromagnetic radiation or Dirac waves describing electrons or neu- 
trinos? Does a perturbed BH have a signature emission like a 'ringing', analogous to the 
case of a struck bell? To answer such questions, Chandra devoted himself to studying BH 
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perturbations from 1970s onwards^^'^*'"^'^. 

When a BH is perturbed, the curved space-time geometry around the BH will be subjected 
to metric fluctuations. For sufficiently small perturbations, a linear analysis of the metric 
ffuctuations can be carried out in terms of normal modes except that dissipation due to 
both emission of gravitational waves as well as their absorption by the BH make the mode 
frequencies complex, with the decay reffected in the imaginary parts. In the case of a 
perturbed BH, such quasi-normal modes (QNMs) correspond to a characteristic ringing of 
the ambient space-time geometry that eventually decays due to emission of gravitational 
waves. 

QNMs were discovered by Vishveshwara^^ and Press^^ while studying gravitational wave 
perturbations of BHs. Chandra and Detweiler suggested for the first time numerical methods 
for calculating the QNM frequencies^^. Such investigations throw light on methods for 
direct detection of BHs. For example, matter falling into a Schwarzschild BH would lead 
to excitation of QNMs, resulting in emission of gravitational waves with a characteristic 
frequency that is inversely proportional to the BH mass. 

One can understand this dependence from simple dimensional analysis. QNMs would in- 
volve perturbations of the event horizon characterized by the Schwarzscild radius Rg (eq.(5)). 
So, the oscillation wavelengths would be typically of a size proportional to Rg, making the 
frequencies depend inversely on the BH mass. A supermassive BH with mass 10^ Mq would 
ring with a frequency of about 10^^ Hz. Because of seismic noise, LIGOs cannot detect 
gravitational waves having such low frequencies. Only a space-based gravitational wave 
detector like LISA (Laser Interferometer Space Antenna) can pick up such low frequency 
signals from supermassive BHs^^'™. 

Chandra developed innovative techniques to study BH perturbations, and showed that 
radial and angular variables could be decoupled to obtain separable solutions for Dirac 
equation in Kerr background, corresponding to a massive particle (like an electron) ^^. Using 
similar techniques, Don Page extended the separation of variables for massive Dirac equation 
to the Kerr-Newman case^^. In 1973, Teukolsky had separated the Dirac equation for two 
component massless neutrinos in the Kerr backgrounds^. It will be interesting to investigate 
if Chandra's technique can succeed in separating the Dirac equation for massive neutrinos 
(with flavour mixing and massive right-handed components included) in the Kerr or Kerr- 
Newman background. 
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Kerr BHs possess ergosphere, a region surrounding the event-horizon where test particles 
with negative angular momenta (i.e. with reverse sense of rotation relative to BH rotation) 
can have negative energy (as measured by a distant inertial observer) orbits. Penrose, in 
1969, had shown an ingenious way to extract rotation energy of a Kerr BH that involved 
sending an object that breaks up into two in the ergosphere, with one of the parts going 
into a negative energy trajectory, while the other escaping with an energy greater than the 
initial energy (since energy is conserved) ^^. 

The wave analogue of Penrose process is superradiance wherein impinging scalar, elec- 
tromagnetic or gravitational waves emerge out with greater energy after scattering off Kerr 
BHs. Zel'dovich was the first to show the existence of superradiance in 1970^^. Chandra 
and Detweiler undertook a thorough investigation of scattering of electromagnetic, gravita- 
tional and neutrino waves in the Kerr background, and showed that neutrinos do not exhibit 
superradiance'^^. Absence of neutrino superradiance is most likely due to PEP^^^^^. 

Exact solutions of two plane gravitational waves colliding with each other were obtained 
for the first time by Szekeres^^ as well as Khan and Penrose^^. Their work showed that due 
to mutual gravitational focusing, the collision leads to curvature singularity where gravity 
becomes infinite. Chandra, along with Valeria Ferrari and Xanthopoulos, showed that the 
mathematical theory of colliding gravitational waves can be cast in the form of mathematical 
theory of BHs, and that the formation of curvature singularity due to gravitational focusing 
IS generic . 

In the later years, Chandra and Valeria Ferrari studied non-radial oscillations of rotating 
stars taking into account general relativistic effects^^"^^. They showed that the oscillations 
could be described in terms of pure metric perturbations, reducing the problem to scattering 
of gravitational waves in curved space-time geometry. For strongly gravitating objects like 
neutron stars, such gravitational waves may get trapped inside due to deep gravitational 
potential well, leading to trapped modes that survive for long durations. 

Chandra was awarded the Nobel prize for physics in 1983. His method of studying 
diverse astrophysical topics involved applying physical theories that had been corroborated 
experimentally, and then subjecting the relevant equations that followed from theory to 
rigorous and innovative mathematical analysis. No wonder that most of the new results he 
obtained were later confirmed by observations. 
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